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LONG-TERM  GOALS 

The  primary  objective  of  the  proposed  research  is  continued  development  of  an  APL-UW  synthetic 
aperture  sonar  simulation  tool  (SAS-ST).  This  tool  can  generate  sets  of  realistic  signals  suitable  for 
SAS  and  acoustic  color  template  (i.e.,  target  strength  as  a  function  of  aspect  angle  and  frequency) 
processing  for  targets  near  the  seafloor.  Updates  to  the  propagation  component  are  tested  and 
validated  against  experimental  data. 

OBJECTIVES 

The  proposed  research  includes  a  few  specific  aims  for  improvements  to  SAS-ST.  First,  a  fast  ray- 
based  model  for  propagation  in  a  homogenous  waveguide,  which  tracks  time-of-flight  wave  packets 
from  sources  to  targets  and  then  to  receivers,  was  reduced  a  fluid-fluid,  half-space  model  and  applied 
to  a  circular  SAS  simulation.  The  ray  model  use  image  sources  and  receivers  to  account  for  interaction 
with  the  water-sediment  interface.  Validation  of  the  half-space  model  was  carried  out  by  comparison 
to  Pond  Experiment  2010  (PondExlO)  data  and  against  finite-element  model  (FEM)  simulations. 

APPROACH 

A  ray-base  model  developed  in  FY 1 1  assumed  a  homogeneous  layer  water  between  an  upper  semi¬ 
infinite  half-space  of  air  and  a  lower  semi-infinite  half-space  of  a  homogenous  sediment.  Rays  are 
assumed  to  travel  in  straight  line  segments  from  a  source  location  to  a  target  location  and  after 
interacting  with  the  target  to  a  receiver  location.  For  an  interaction  of  a  ray  with  an  interface,  an 
appropriate  reflection  coefficient  is  included  in  the  contribution  to  the  scattered  signal.  At  present  the 
sediment  can  be  modeled  as  either  an  attenuating  fluid  with  a  frequency-independent  loss  parameter  or 
a  effective  density  fluid  model.  For  the  rays  that  interact  with  the  interfaces,  these  rays  can  be 
represented  as  coming  from  image  sources  or  arriving  at  image  receivers.  For  a  (image)  source  or 
(image)  receiver  located  at  r,  and  a  target  located  at  rt,  the  horizontal  and  total  separation  distances  are 
\Rj  -  R,\  and  \r,  -  rt\,  respectively.  The  contribution  of  the  z'th  source  and  the  /th  receiver  to  the 
frequency  spectrum  of  the  total  scattered  signal  can  be  written  as 
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where  Psrc((o)  is  the  frequency  spectrum  of  the  transmitted  wave  packet.  The  reflection  coefficients  at 
the  upper  and  lower  interfaces  are  A(0,j  and  B(Q,).  Here,  8,  simply  denotes  a  local  grazing  angle  at  the 
interface,  which  is  given  by  cos(0;)  =  | R,  -  R\  /  \rt  -  rt |,  and  the  in  and  n  exponents  indicate  the  number 
of  interactions  a  ray  has  with  a  given  interface.  The  time  delay  for  propagation  from  the  source  to  the 
target  is  ti  =  \rt  -  r,|  and  the  delay  from  the  target  to  the  receiver  is  tj  =  | r,-  -  rt |.  The  scattering  process  is 
a  convolution  of  a  free-field  scattering  fonn  function  F(0y,(j)y,co)  with  the  incident  pressure  spectrum 
Psrc( ro).  Here,  0(/  and  <j);/-  are  target-centered  polar  and  azimuthal  angles,  which  are  related  to  the  0,. 

WORK  COMPLETED 

During  FY12,  the  waveguide  model  given  in  (1)  was  reduced  to  a  half-space  model  for  co-located 
source  and  receiver.  This  model  is  suitable  for  monostatic  SAS  application  and  consists  of  the  source 
and  receiver  locations  rs  and  rr  and  a  single  image  source  located  at  rsi  and  image  receiver  at  rri.  This 
yields  four  possible  rays  paths:  S— >T— >R,  S— >B— >T— >R,  S— »T— »B— >R,  and  S— »B— »T— »B— >R,  where 
the  symbol  >”  represents  a  propagation  phase  delay  with  S(ource),  T(arget),  B(ottom),  and 
R(eceiver).  The  path  segments  S— »B— >T  and  T^B— +R  represent  propagation  from  an  image  source  to 
a  target  and  from  a  target  to  an  image  receiver.  It  is  noted  that  the  first  path  is  a  direct  path;  the  second 
and  third  paths  interact  with  the  bottom  once,  scatter  is  in  a  bistatic  direction,  and  are  reciprocal  paths; 
and  the  last  path  is  a  backscattering  path  with  two  bottom  interactions.  Using  d\  =  | rs  -  r, |,  di  =  | rt  -  r,  |, 
d-i  =  | rsi  -  rt |,  =  | r,  -  rri |,  Eq.  (1)  reduces  to 
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with  t\  =  {d i  +  di)  /  c,  t2  =  (di  +  d 3)  /  c,h  =  {d\  +  J4)  /  c,  and  =  (di  +  d$)  /  c;  c  being  the  speed  of 
sound  in  water.  The  scattering  form  functions  -F*(0y, <(>/,•,  co)  in  Eq.  (2)  depend  on  the  locations  of  the 
sources,  receivers,  and  target.  An  inverse  Fouier  transfonn  of  R(co)  then  gives  a  generated  sonar  signal 
that  includes  the  four  primary  acoustics  paths  for  a  target  near  an  interface. 

RESULTS 

The  half-space  propagation  model  was  applied  to  the  problem  of  circular  SAS  with  a  target  located 
within  the  CSAS  path.  For  the  results  presented  here,  the  target  is  placed  at  the  center  of  the  circular 
path.  The  target  is  in  a  broad-side  orientation  at  0°,  180°,  and  360°.  The  tail  of  the  replica  is  pointing 
at  the  source/receiver  at  90°  and  the  nose  is  pointing  at  the  source/receiver  at  270°.  The  radius  of  the 
path  was  varied  from  5  to  50  m  in  0.5  m  increments,  and  adjacent  signals  were  recorded  at  1° 
increments  along  the  circular  path.  The  source  and  receiver  were  located  3.8  m  above  the  sediment 
and  were  assumed  to  be  omni-directional.  The  simulations  consider  the  scattering  from  an  aluminum 
replica  of  a  real  bullet-shaped  artillery  shell.  The  replica  is  ~4  inches  in  diamater  and  ~16  inches  in 
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length  with  its  axis  of  symmetry  parallel  to  the  water-sediment  interface.  This  target  was  deployed 
during  PondExlO  and  finite-element  model  results  are  available  for  comparisons.  The  sound  speed 
and  density  of  the  water  are  c  =  1500  m/s  and  p  =  1000  kg/m  .  For  the  simulations  shown  in  the 
figures  below,  the  sediment  was  considered  to  be  an  attenuating  fluid.  The  sediment  properties  are  cs  = 

■5 

1775  m/s,  p.s  =  2000  kg/m  ,  and  5  =  0.008  (i.e.,  the  loss  parameter).  For  the  simulations,  the  source 
emitted  a  6  ms  chirp  with  a  carrier  frequency  of  16  kHz  and  bandwidth  of  30  kHz.  These  parameters 
roughly  correspond  to  experimental  conditions  during  the  deployment  of  the  APF-UW’s  rail/tower 
during  PondExlO.  Figure  1  shows  the  acoustics  color  templates  obtained  from  PondExlO  data,  FEM 
simulations,  and  the  ray-based,  half-space  model  in  SAS-ST.  In  Eq.  (2),  we  use  a  free-field  scattering 
form  function  derived  from  the  FEM  prediction  of  the  target  response  where  a  Kirchoff-Helmholtz 
integral  is  used  to  propagate  the  scattered  pressure  to  a  10  m  range  from  the  target.  This  scattered 
pressure  is  converted  to  a  scattering  form  function  and  placed  in  a  table:  F(6i,(j)j,c0k).  Given  the  source 
and  receiver  locations  along  the  circular  path  and  the  local  grazing  for  the  target  location,  values  of  9 
and  cj)  are  determined.  Linear  interpolation  was  then  employed  to  extract  scattering  fonn  function 
values  from  the  table.  Inspection  of  Figure  1  shows  that  both  the  FEM  and  SAS-ST  results  capture 
much  of  the  structure  observed  in  the  PondExlO  data.  This  structure  includes  both  geometrically 
reflected  acoustic  energy  and  acoustic  energy  that  couples  into  the  elastic  response  of  the  target. 
Comparison  of  Figure  1(b)  and  1(c),  particularly  in  the  -10  to  -17  dB  range  reveals  that  the  SAS-ST 
result  has  a  slight  smoothing.  This  is  a  consequence  of  the  discrete  representation  of  ^(0i,4>j,cok)  and 
the  linear  interpolation. 
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Figure  1:  Acoustic  color  templates  for  an  aluminum  replica  of  a  bullet-shaped  artillery  shell  at  a  10 
m  range,  (a)  PondExlO  measurement,  (b)  FEM  and  Kirchoff-Helmholtz  integral  result,  (c)  SAS- 
ST’s  ray-based,  half-space  model  using  a  scattering  form  function  derived  from  the 

scattered  pressure  computed  in  (b). 
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Figure  2:  Acoustic  color  templates  for  an  aluminum  replica  of  a  bullet-shaped  artillery  shell  at  a  5 
m  range,  (a)  PondExlO  data,  (b)  FEM  model  and  Kirchoff-Helmholtz  integral  result,  (c)  SAS- 

ST ’s  ray-based,  half-space  model. 


Figure  2  shows  the  PondExlO  data  for  the  aluminum  replica,  FEM  and  Kirchoff-FIelmholtz  integral 
simulation,  and  the  SAS-ST’s  ray-based,  half-space  model  for  the  target  at  a  5  m  range.  An  important 
aspect  of  the  SAS-ST  simulation  is  that  it  used  a  F(9i,(j)j,cok)  that  was  derived  from  a  free-field  FEM 
and  Kirchoff-Helmholtz  integral  simulation  for  a  target  at  a  10  m  range.  This  is  possible  because  the 
scattering  form  function  is  range  independent  and  contains  only  the  angular  infonnation  and  the 
construction  of  the  model  in  Eq.  (2)  explicitly  accounts  for  interactions  with  the  water-sediment 
interface.  Thus,  only  a  single  free-field  measurement  from  a  target  or  a  single  free-field  FEM 
simulation  is  sufficent  to  provide  a  scattering  form  function  that  can  be  used  to  rapidly  predict  the 
scattering  from  a  target  near  a  bound  at  any  range  (provide  the  range  meets  the  far  field  requirement). 

IMPACT/APPLICATIONS 

The  current  version  of  SAS-ST  has  been  used  in  comparisons  with  data  obtained  during  PondexlO  and 
Pond  Experiment  2009.  In  addition,  the  new  half-space  model  included  in  SAS-ST  a  significantly 
reduced  computational  burden  (i.e.,  both  memory  and  execution  time)  over  high-fidelity  FEM  and 
Kirchoff-Helmholtz  integral  methods.  The  new  half-space  model  permits  one  to  investigate 
operational  patterns  in  MCM  scenarios  such  as  circular  SAS  and  can  be  used  to  generate  sets  of  SAS 
signals  in  testing  automatic  target  recoginition  algorithm. 

RELATED  PROJECTS 

The  companion  award,  “Generation  of  Synthetic  SAS  Data  for  Targets  near  the  Seafloor:  Target 
Scattering  Component,”  seeks  to  improve  scattering  models  for  targets  near  an  interface  used  in 
SASSSY.  The  research  described  in  the  present  Annual  Report  will  be  merged  with  the  target 
scattering  component  to  provide  an  overall  improvement  to  the  fidelity  of  the  simulated  SAS  data. 
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